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We study the cross-correlations in stock price changes between the S&P 500 companies by 
introducing a weighted random graph, where aU vertices (companies) are fuUy connected, and each 
edge is weighted. The weight assigned to each edge is given by the normalized covariance of the 
two modified returns connected, so that it is ranged from -1 to 1. Here the modified return means 
the deviation of a return from its average over all companies. We define infiuence-strength at each 
vertex as the sum of the weights on the edges incident upon that vertex. Then we found that the 
infiuence-strength distribution in its absolute magnitude \s\ follows a power-law, -P(|s|) ~ |s|~ , with 
exponent 5 ~ 1.8(1). 

PACS numbers: 89.75.Da, 89.65.Gh, 05.10.-a 
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Recently complex systems such as biological, economic, 
physical, and social systems have received considerable 
attentions as an interdisciplinary subject Q. Such sys- 
tems consist of many constituents such as individuals, 
companies, substrates, spins, etc, exhibiting cooperative 
and adaptive phenomena through diverse interactions 
between them. In particular, in economic systems, adap- 
tive behaviors of individuals, companies, or nations, 
play a crucial role in forming macroscopic patterns such 
as commodity prices, stock prices, exchange rates, etc, 
which are formed mostly in a self-organized way B . Re- 
cently, many attentions and studies have been focused 
and performed on the fluctuations and the correlations 
in stock price changes between different companies in 
physics communities by applying physics concepts and 
methods |,|. 

Stock price changes of individual companies are in- 
fluenced by others. Thus, one of the most important 
quantities in understanding the cooperative behavior in 
stock market is the cross-correlation coefficients between 
different companies. Since the stock prices changes de- 
pend on various economic environments, it is very hard 
to construct a dynamic equation, and predict the evolu- 
tion of the stock price change in the future. Recently, 
there have been some efforts to understand the correla- 
tions in stock price changes between different companies 
using a random matrix theory, where large eigenvalues 
are located far away from the rest part predicted by the 
random matrix theory, reflecting the collective behavior 
of the entire market as well as the subordination among 
the stock prices 



Let Yi{t) be the stock-price of a company i (i = 
1, . . . , N) at time t. Then the return of the stock-price 
after a time interval At is defined as 

S,it)=lnY,{t + At)-\nY,it), (1) 

meaning the geometrical change of Yi(t) during the in- 



terval At. We take At = 1 day for the following analysis 
in this Letter. The cross-correlations between individ- 
ual stocks are considered in terms of a matrix C, whose 
elements are given as 
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where the brackets mean a temporal average over the 
period we studied. Then C^j can vary between [-1,1], 
where Cij = 1 (-1) means that two companies i and j are 



completely correlated (anti-correlated), while Ci^ 







means that they are uncorrelated. Since the matrix C 
is symmetric and real, all eigenvalues are real, and the 
largest eigenvalue is not degenerate. It was found that 
the eigenvector corresponding to the largest eigenvalue 
is strongly localized at a few companies which strongly 
influence other companies in the stock price changes 
[H,g. Moreover, the ultrametric hierarchical tree struc- 
ture was constructed among those companies, using the 
concept of the minimum spanning tree in the graph the- 
ory, which arranges those companies in order following 
their strengths of influence MM . 

In this Letter, we study further properties of the cross- 
correlations in stock price changes using a random graph 
theory |10). The study of complex systems through a 
random graph was initiated by Erdos and Renyi (ER) 
[ pdl , however, the ER model is too random to describe 
real-world complex systems. Recently, Barabasi and Al- 
bert (BA) [|l2[ introduced an evolving network where the 
number of vertices increases with time rather than fixed, 
and a newly introduced vertex is connected to m already 
existing vertices, following the so-called preferential at- 
tachment rule. The number of edges k incident upon a 
vertex is called the degree of the vertex. Then the degree 
distribution P{k) follows a power-law P{k) ^ k~'^ with 
7 = 3 for the BA model, while for the ER model, it 



follows a Poisson distribution. The network following a 
power-law in the degree distribution is called scale-free 
(SF) network. 

For the problem of the correlations in stock price 
changes, each vertex (edge) in the random graph corre- 
sponds to a company (the cross-correlation in stock price 
changes between the companies connected via that edge). 
The random graph generated in this way is different from 
a typical one in the following way: While the edge in a 
typical random graph has weight either 1 or 0, depend- 
ing on if the edge is present or not, respectively, the 
edge in the random graph we will introduce has weight 
{wij}, which is rather distributed between [-1,1]. For 
further studies, the random graph for the former (later) 
case is called binary random graph (BRG) (weighted 
random graph (WRG)). While WRG can be found easily 
in real- world networks such as neural networks, cardio- 
vascular networks, and respiratory networks in biological 
systems, acquaintance networks in social g^^tems, etc, it 
is less studied, compared with BRG ||l^,|lj|. In WRG, 
one may wonder if the edge with weak weight, called 
the weak edge, can be ignored in considering the cor- 
relations, however, there have been ongoing discussions 
about the importance of the weak edges, for example, in 
social acquaintance network, the scientific collaboration 
web, and ecosystems [n5]jlq|. We will also show that 
the WRG introduced in this Letter is different from the 
corresponding DRG in which weak edges are ignored and 
the remaining edges are assigned a unit weight. 

Recently, Yook, Jeong and Barabasi (YJB) introduced 
a WRG p6[ . In this model, a vertex i is newly introduced 
at each time step, connecting to m vertices existing al- 
ready according to the so-called preferential attachment 
rule. The edge connecting from the vertex i to an exist- 
ing vertex j is assigned a weight Wij, depending on the 
degree of the vertex j, which is not uniform. The weight 
at each vertex is defined as the sum of the weights on the 
edges incident upon that vertex, which follows a power- 
law in its distribution, P{w) ~ w^ , where w means 
weight at a vertex. The exponent 6 is different from the 
degree exponent 7, and turns out to strongly depend on 
the mean degree m. While the YJB model is meaningful 
as a first step towards understanding diverse WRGs in 
real- world in a simple way, it is very different from what 
we will consider in association with the cross-correlations 
in stock price changes. The WRG we will introduce be- 
low are fully connected, and each edge is weighted, while 
the YJB graph is also weighted, but sparsely connected. 

We consider the cross-correlations in stock-price 
changes between the S&P 500 companies during 5-year 
period 1993-1997. Thus, the N = 500 companies cor- 
respond to 500 vertices, which are fully connected to 
each other through N(N — l)/2 edges. Each edge is as- 
signed a weight, Wij {i,j = 1, . . . , N), which is slightly 
modified from the cross-correlation coefficient Cij de- 



fined in Eq.(2). Before defining Wi^j specifically, we first 
recall some properties exhibited by Cij. It is known 
that the distribution of the coefficients {Ci,j} is a bell- 
shaped curve, and the mean value of the distribution is 
slowly time-dependent, while the standard deviation is 
almost constant [0 . The time-dependence of the mean 
value might be caused by external economic environ- 
ments such as bank interest, infiation index, exchange 
rate, etc, which fluctuates from time to time. To ex- 
tract intrinsic properties of the correlations in stock price 
changes, we introduce a quantity. 



G,(t) = S,(t)-l5]^,(t) 



(3) 



where Gi (t) indicates the relative return of a company i 
to its mean value over the entire 500 companies at time t. 
The cross-correlation coefficients are redefined in terms 
oiG,, 
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The cross-correlation coefficients {wij} are assigned to 
each edge of the WRG as its weight. In order to check 
if the distribution P{w) is time- independent, we take 
temporal average in Eq.(4) over each year from 1993 to 
1997. In Fig.|l|, we plot the distributions of {wij} ob- 
tained for each year. The data for different years are 
indeed overlapped, and time-independent. Therefore we 
think that the cross-correlation coefficients {wij} we 
introduced are appropriate to study intrinsic properties 
of the cross-correlations between the 500 S&P companies. 

We define the influence-strength Si at a vertex i as the 
sum of the weights on the edges incident upon the vertex 
i, that is. 
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where j denotes the vertices connected to the vertex i. 
Here {wij} was obtained numerically by temporal av- 
eraging over the 5 years in Eq.(4). Then the weight Si 
means the net amount of influence for the company i to 
affect other companies in stock price changes. Since the 
weight Wij is distributed between [-1,1], the inffuence- 
strength at a certain vertex could be negative. Thus, 
we deal with the absolute magnitude of the influence- 
strength for each vertex. In Figjl, we plot the influence- 
strength distribution -P(|s|) in the absolute magnitude 
as a function of |s|, which turns out to follow a power- 
law, -P(|s|) ~ \s\~^ ■ The exponent S is estimated to 
be (5 « 1.8(1). Thus the cross-correlations in stock 
price changes forms a SF network, which in particular 
is a weighted SF network. To our knowledge, this is 
the first observation of SF WRG emerging in real-world 
economic systems. The presence of a scaling in the 



cross-correlations through the WRG would be related 
to the ultrametric hierarchical tree structure, implying 
that there exist a few companies having strong influence 
in stock price changes. On the other hand, it is easy 
to know that as the degree exponent in SF networks 
is smaller, the connectivity to the hub, the vertex with 
the largest degree, is higher, and the network is much 
centralized. This fact is also applicable to the SF WRG 
we introduced. Since the influence-strength exponent 6 
is smaller than 2 in the WRG, the vertex having largest 
influence-strength plays a much important role in af- 
fecting stock price changes of other vertices, compared 
with the role of the hub in the Internet [|l8| , the world- 
wide web [Q and the metabolic network pCI| ] , where the 
degree exponent is greater than 2. We think that this 
result reflects that economic systems is much correlated 
and adaptive to achieve high proflts. Thus the rich-get- 
richer phenomenon appears much strongly in economic 
systems, compared with the information systems. In 
contrast, we can expect that a simple drop in the stock 
price occurring in one of most influential companies could 
lead to a crash in the entire stock market. 

Next, we study the matrix W whose elements are given 
by Wi.j. We found the distribution of the eigenvalue A 
for the matrix W shown in Fig.^, lookingsimilar to the 
one for the matrix C previously studied MM; For small 
A, the spectrum is likely to follow the prediction by the 
random matrix theory, while it exhibits a fat-tail behav- 
ior for large A. However, we found that the gap between 
the first and second largest eigenvalues for the matrix 
W is not as big as the one for the matrix C. Thus we 
think that the huge gap observed in the eigenspectrum 
for the matrix C may be caused by some external effects. 

Finally, we also consider the components {wj.i} of the 
eigenvector corresponding to the largest eigenvalue for 
the matrix W, where the vertices are ordered following 
the absolute magnitude of influence-strength. Accord- 
ing to the matrix theory, the square of the component, 
{v^ i}, indicates the relative strength of the contributions 
for each vertex to the largest eigenvalue. We plot the 
square of the component of the eigenvector {v^ i } versus 
the index j in Fig. 4. It is found that the components 
{v^ i} is strongly localized at the vertices with strong 
influence-strength as shown in Fig. 4. Thus, the compa- 
nies with large influence-strength affect other companies 
strongly in stock price changes. Interestingly, we found 
another SF behavior, w^ ^ ~ j~^ with (3 « 1-8(2) for large 
J, as shown in the inset of Fig. 4, which is consistent with 
the one S « 1.8(1), obtained from the influence-strength 
distribution. 

It would be interesting to compare our result with 
the one by Vandewalle et al. |2l|, who eliminated weak 
edges, and regarded the remaining sparse random graph 
as a BRG. To be specific, they considered the cross- 
correlation coefficients defined in Eq.(2) in stock price 



changes between the 6 358 US companies, and assigned 
them to each edge of the fully connected graph. Next, 
they chose N strongest edges out of N{N — l)/2 edges 
in the magnitude of Cij, and assigned a unit strength 
to each of them, so that weak edges are ignored. Then 
the graph becomes sparse and binary. They measured 
the degree distribution for this BRG, which follows a 
power-law, P{k) '^ k~^ ^ with 7 « 2.2, which is obviously 
different from our results, 5 ~ 1.8(1) and (3 « 1.8(2). 



In conclusions, we have considered the cross- 
correlations in stock price changes between the S&P 
500 companies by introducing a weighted random graph 
(WRG). The vertices of the WRG representing the 500 
companies are fully connected to each other through 
weighted edges. An edge connecting vertices i and j has 
the weight given by the correlation coefficient Wij, de- 
fined as the normalized covariance of modified returns of 
two companies i and j. Here the modified return of a 
company i means the deviation of the return of the com- 
pany i from its average over all 500 companies. This mod- 
ification yields the effect of excluding the overall behavior 
of the entire stock prices fluctuating from time to time. 
The distribution of the correlation coefficients obtained 
using the modified return is time- independent, and the 
coefficients themselves describe generic correlations be- 
tween different companies without considering the effect 
of external environments. Next, we defined the influence- 
strength at each vertex as the sum of the weights assigned 
to the edges incident upon that vertex. We found that 
the influence-strength distribution follows a power-law 
P(|s|) ^ |s|"'' with 5 « 1.8(1), where s means influence- 
strength. The fact that the exponent 5 is smaller than 2 
implies that the stock price changes of the 500 companies 
are much strongly correlated, compared with the Internet 
topology, or the world-wide web, reflecting that cooper- 
ative and adaptive phenomena appear much dominantly 
in economic systems. 
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FIG. 2. Plot of the influence-strength distribution -P(|s[) 
versus the absolute magnitude of the influence-strength |s| 
in double-logarithmic scales. The solid guideline has a slope 
-1.8. 
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FIG. 3. Plot of the distribution of the eigenvalues of the 
matrix W. Inset: Plot of the distribution of the eigenvalues 
versus eigenvalue in the fat-tail region. 



FIG. 1. Plot of the distribution of the cross-correlation co- 
efficients {wij}. The data are obtained by temporal averaging 
over each year from 1993 to 1997. 
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FIG. 4. Plot of the component squares v'ji of the eigen- 
function corresponding to the largest eigenvalue versus the 
vertex index j. The vertices are ordered following the abso- 
lute magnitudes of their influence-strengths. Inset: The same 
plot as the main panel but in double-logarithmic scales using 
log-bin. The straigth guideline represents v^ ^ ~ j^^'^- 
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